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16  Abstract 

Criteria  are  presented  for  [selection  of  analytic  models  to 
represent  waiting  times  due  to  queing  processes . An  existing 
computer  model  by  M.F.  Neuts  which  assumes  general  nonparametr ic 
distributions  of  arrivals  per  unit  time  and  service  times  for  a 
single  service  was  envisioned  as  best  fulfilling  requirements. 

Data  obtained  from  Denver  Stapleton  Airport  were  applied  to  this 
model.  Service  ti.rn.es  and  arrival  rates  at  an  express  baggage  check 
facility,  a security  station,  and  a gate  were  used  as  inputs.  Delay 
times  corresponding  to  the  observed  arrival  rates  were  recorded  and 
compared  to  model  outputs.  Using  the  T-test,  agreement  was  obtained 
at  the  5 percent  level  of  significance  for  the  mean  values  of  the  first 
two  facilities.  Predictions  of  waiting  time  distribution,  however, 
did  not  pass  the  Kolmogroff  test  at  the  same  level  of  significance. 
Discrepancies  are  due  to  a lack  of  time  resolution  in  arrival  times 
and  the  application  of  this  model  to  multiserver  situations. 
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PREFACE 


This  study  was  undertaken  as  part  of  Project  Plan  Agreement 
FA-632.  One  objective  of  this  program  is  the  development  and  test- 
ing of  techniques  to  evaluate  airport  landside  congestion  param- 
eters such  as  delay.  This  analysis  provides  a method  of  determin- 
ing delays  arising  from  queuing.  Results  contained  in  this  report 
should  be  considered  as  preliminary  and  indicative  of  further 
development,  especially  where  multi-server  landside  facilities 
require  analysis. 

The  authors  would  like  to  express  their  appreciation  to  the 
Program  Manager,  Mr.  Murk  Gorstein,  for  his  support  and  con- 
structive criticsm. 
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EXECUTIVE  SUMMARY 

An  analytical  queuing  model  was  examined  for  accuracy  in 
predicting  one  of  the  major  measures  of  airport  congestion,  namely, 
the  waiting  time  at  a landside  facility.  The  particular  model 
chosen  is  applicable  to  a single  channel  server  whose  arrival  and 
service  time  processes  have  general  distributions.  This  model  was 
selected  because 

a.  Arrival  time  distributions  obtained  from  measurements  were 
tested  and  found  to  be  non-Poisson. 

b.  A potentially  wide  range  of  applicability  was  desired 
because  of  the  many  service  facilities  present  at  the  landside. 

c.  This  model  was  available  as  a computer  program  which 
furnished  waiting  time  and  queue  length  distributions  directly. 

In  order  to  investigate  the  feasibility  of  using  this  queuing 
model  as  a tool  for  evaluating  landside  congestion,  statistical 
comparisons  were  performed  to  determine  whether  there  is  any 
correspondence  between  observed  field  data  and  model  predictions. 
Using  observed  passenger  arrival  rates  and  service  time  distribu- 
tions obtained  at  Denver-Stapleton  Airport  as  the  inputs,  the  model 
validation  was  done  by  comparing  predicted  waiting  times  at  the 
security  station,  express  check-in,  and  one  boarding  gate.  Using 
the  t-test,  agreement  was  obtained  at  the  5 percent  level  of  sig- 
nificance for  the  mean  values  of  the  first  two  facilities.  More 
comprehensive  data  collection  is  required  to  validate  the  probabil- 
ity distributions  of  the  waiting  time  and  queue  length. 

This  single  channel  queuing  model  may  be  viewed  as  the  first 
step  in  evaluating  landside  congestion  using  an  analytical  approach. 
Further  developments  in  this  area  would  include  multi-channel  facil- 
ity modeling  and  the  linking  of  various  facilities  to  formulate  a 
congestion  model  for  the  complete  landside. 
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1.  INTRODUCTION 


Some  of  the  Denver  - Stapleton  International  Airport  passenger 
data  which  were  collected  for  simulation  model  validation  were 
used  to  determine  the  validity  of  an  analytic  model.  This  was 
done  by  comparing  the  analytically  computed  queue  length  and  wait- 
ing time  distributions  with  the  actual  observed  queue  length  and 
waiting  time  distributions  in  the  field.  It  was  hoped  that  this 
validation  work  would  lead  to  the  determination  of  the  feasibility 
of  using  analytic  models  for  estimating  the  level -of -service  of  the 
landside  facility. 

The  input  data  for  the  analytic  model  are  the  actual  arrival 
rates  and  the  actual  service  time  as  observed  in  Stapleton  Inter- 
national Airport  in  Denver.  The  outputs  of  the  model  are  the  queue 
length  and  waiting  time  distributions  computed  for  each  facility 
under  study.  The  waiting  time  output  was  compared  with  its  field 
observed  counterpart  for  model  validation.  It  is  believed  that 
this  quantity  is  one  of  the  most  meaningful  indicators  of  the  level 
of  service. 
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2.  THE  QUEUING  MODEL 


2.1  MODEL  SELECTION 

There  are  six  basic  characteristics  which  can  be  used  to 
specify  a queuing  model: 

1.  Arrival  Pattern:  The  arrival  pattern  of  the  "customers" 

could  be  deterministic  (D)  or  described  parametrically  as  a Poisson 
Distribution  (or  exponential  arrival  time)  (M)  or  general  (0). 

For  many  airport  landside  facilities  the  arrival  pattern  is  most 
likely  general  (G) . 

2.  Service  Pattern:  The  service  time  pattern  of  each  server 

channel  could  also  be  deterministic  (D) , exponential  (M) , or 
general  (G) . In  the  airport  landside  situation  it  is  most  likely 
general  (G) . 

3.  Number  of  Service  Channels:  The  queuing  process  may  take 

place  at  a facility  with  a single  channel  or  a number  of  parallel 
channels.  In  the  airport  landside  both  the  single  and  parallel 
channel  cases  are  applicable. 

4.  System  Capacity:  Facilities  may  accommodate  a limited 

queue  size  or  the  maximum  queue  size  may  be  unbounded.  Since  it 
is  most  likely  that  no  airport  landside  customers  will  be  turned 
back,  the  maximum  possible  queue  length  may  be  considered  infinite 
for  most  queuing  situations  of  interest. 

5.  Number  of  Service  Stages  (in  series  per  each  channel): 

Each  channel  may  have  a single  stage  or  multiple  stages  connected 
in  series.  Both  single  stage  and  multi-stage  processes  are 
manifested  by  the  airport  landside. 

6.  Queue  Discipline:  Queue  discipline  can  be  classified  into 

(a)  f i r s t -c ome - f i r s t - s erved  (ICES),  alias  f i rst - i n- f i rst -out 
(FIFO);  (b)  service  in  random  (SIRO);  (c)  1 as t-come-f i rst-served 
(LCFS)  ; Ed)  priority  (PR1);  and  (e)  general  discipline  (GD)  . The 
airport  landside  system  operates  mostly  on  a f i rst-comc- f i r st - 
served  (FCFS  or  FIFO)  basis,  with  some  priority  (PRI)  treatments. 


Based  on  the  six  basic  characteristics,  the  most  desirable 
model  for  the  airport  landside  system  would  be  a model  with 
(1)  general  arrival  pattern,  (2)  general  service  time  pattern, 

(3)  multiple  channels,  (4)  infinite  or  finite  maximum  queue  length, 
(5)  single  or  multiple  stages,  per  channel,  and  (6)  first-come- 
first-served  and/or  priority  service  discipline.  In  queuing  theory 
notation  this  most  desirable  model  is  represented  as: 


G/G/Nj/ 


with  multiple-stage  per  channel 


It  should  be  realized  that  the  so  called  most  desirable  model  is 
probably  the  most  complicated  model  one  would  ever  attempt  to 
utilize. 

The  development  of  analytic  queuing  models  appear  to  have  been 
influenced  strongly  by  two  factors:  the  origins  of  models  in  the 
study  of  congestion  in  telephone  systems  and  the  question  of  what  is 
easy  and  possible  in  mathematical  analysis.  Therefore,  many  assump- 
tions commonly  made  in  queuing  analyses  are  precisely  those  which 
seem  reasonable  in  constructing  queuing  models  associated  with  tele- 
phone systems  and  those  which  make  the  models  mathematically  solvable. 
As  a result,  analytic  models  which  deal  with  multiple  channels  or 
servers  are  generally  restricted  to  Poisson  arrivals  and  exponential 
service  times  fi.e.,  M/M/N) . Models  which  deal  with  general 
arrivals  and  general  service  times  are  presently  limited  to  single 
channel  ( < 1 / G / 1 ) cases. 

The  simplifying  assumptions  used  in  M/M/N  model  development 
are  not  always  applicable  to  queues  other  than  those  encountered  in 
telephone  systems;  they  are  particularly  inappropriate  when  the 
customers  and  servers  are  really  people  as  in  the  airport  landside. 

In  a telephone  system  the  process  of  switching  channels  or  jockey- 
ing is  carried  out  automatically  and  rapidly,  whereas  in  many  multi- 
counter operations  there  are  no  easily  defined  rules,  and  people 
move  with  inertia.  This  was  observed  visually  at  Denver -Stapleton 
International  Airport  where  queues  were  formed  at  ticket  counters 
at  the  same  time  that  servers  were  observed  idle.  This  point  was 
further  verified  numerically  by  the  field  data.  Out  of  five  sets 
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of  arrival  data  tested,  only  one  set  was  Poisson,  and  none  of  the 
service  times  observed  were  exponential.  This  was  the  main  reason 
why  G/C’./l  model  was  selected  over  M/M/ N for  validation.  Although 
this  single  channel  model  can  represent  only  a limited  number  of 
facilities  in  the  airport  landside,  it  was  chosen  to  determine 
whether  there  is  any  resemblance  between  observed  field  data  and 
the  model  predictions.  In  addition,  a multiple -channel  G/G/N 
model  is  currently  under  investigation. 

2.2  NF.UTS  G/G/l  MODEL 

A readily  available  G/G/l  queuing  model  was  chosen  for  model 
validation.  This  model,  developed  by  Professor  Marcel  F.  Neuts 
of  Purdue  University,  accepts  general  arrivals  and  general  service 
times  and  computes  the  time  dependency  features  of  a single  server 
discrete  time  queue  with  a finite  queue  length.  The  model  was 
coded  in  Fortran  IV  and  it  resides  on  the  DEC-10  computer  system 
at  Transportation  Systems  Center/U.S.  Department  of  Transportation 
(TSG) , Cambridge,  Mass.  The  program  listing  is  included  in  the 
append  ix . 

In  the  model,  the  numbers  of  arrivals  during  successive  unit 
time  intervals  are  independent,  identically  distributed  random 
variables.  These  are  expressed  as  the  probability  of  a number  of 
customers  joining  the  queue  during  a specified  unit  time.  The 
service  times  of  successive  customers  are  also  independent  and 
identically  distributed  random  variables  with  a given  probability 
density.  The  queue  discipline  is  ins. 

For  any  given  time,  n,  the  queue  length  is  denoted  as  X^  and 

the  residual  service  time  (i.e.,  the  number  of  additional  units  of 

service  time  required  by  the  customer  in  service)  is  denoted  by 

Yn . The  joint  probability  density  function  of  queue  length  Xn  and 

residual  service  time  Y„  is  denoted  as  P fi,j),  where 

n n J 

■ p {*„  ■ >■  - j}  • 

Recurrence  relations  were  derived  which  compute  Pn+^(i,j)  from 
P (i,j)  under  all  possible  conditions  of  i and  j.  This  joint 
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probability  density  function  is  used  to  compute  the  distributions 
of  queue  length  and  waiting  time.  These  are  the  basic  outputs  of 
the  model.*  As  stated  before,  the  model  has  been  translated  into 
a Fortran  IV  program.  The  inputs  to  the  program  are: 

a.  The  size  of  the  waiting  room  LI.  In  the  program,  LI  is 
limited  to  100  customers  in  the  system. 

b.  The  maximum  number,  K,  of  arrivals  per  unit  time.  In  the 
program,  K is  limited  to  be  less  than  LI  arrivals. 

c.  The  probability,  P(J),  that  J customers  join  the  queue 
per  unit  time,  where  J = 1,2,  ...  K. 

d.  The  maximum  service  time,  L2,  in  terms  of  time  units.  In 
the  program,  L2  is  limited  to  be  less  than  31  time  units. 

e.  The  probability,  R(J),  that  a service  time  lasts  for  J 
units  of  time,  where  J = 1,  2,  ...  L2. 

f.  The  probability,  PO , that  no  customers  arrive  during  a unit 
of  time. 

g.  The  initial  queue  length,  10,  and  the  initial  residual 
service  time,  JO. 

h.  The  program  control  parameters,  including  the  number  of 
time  units  to  be  run  and  the  output  option  selection  parameters. 

The  full  output  of  this  program  includes  the  following: 

a.  The  mean  queue  length  as  a function  of  time. 

b.  The  distribution  of  queue  length  as  a function  of  time. 

c.  The  mean  waiting  time  as  a function  of  time. 

d.  The  distribution  of  waiting  time  as  a function  of  time. 

e.  The  joint  density  of  the  queue  length  and  the  residual 
service  time  as  a function  of  time. 


For  details,  see  M.F.  Nuets,  "The  Single  Server  Queue  in  Discrete 
Time -Numerical  Analysis,  I,"  Naval  Research  Logistics  Quarterly, 

20  (1973),  297-304. 
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Any  one  or  the  combination  of  the  above  output  items  can  be  printed 
out  according  to  the  user-selected  options. 


2.3 


SUPPORTING  PROGRAMS 


In  addition  to  the  queuing  model  itself,  two  supporting  pro- 
grams were  written.  One  of  them  is  used  to  process  the  raw  data 
for  queuing  model  input.  The  other  program  is  used  to  compare  the 
model-predicted  waiting  time  distribution  with  the  field-observed 
waiting  time  distribution. 


2.3.1  Preprocessing  Program  CH I SQ . F4 

This  program  was  written  for  preprocessing  the  raw  arrival 
data  which  were  gathered  in  Stapleton  Airport.  The  functions  of 
this  program  are: 

a.  Convert  observed  arrivals  into  an  arrival  distribution 
( h i s tog  ram) . 

b.  Determine  the  sample  mean  and  variance  of  the  data. 

c.  Perform  X"-test  (chi  square  test)  on  raw  data  to  determine 
whether  the  arrivals  are  from  a Poisson  process. 

d.  Convert  the  arrival  distribution  into  compatible  units  for 
the  G/G/ 1 mode  1 . 

This  program  was  necessary  because  the  raw  data  were  not  directly 
applicable  to  the  G/G/l  model.  The  test  of  Poisson  was  incorpor- 
ated in  this  program,  because  if  the  majority  of  the  arrivals  were 
Poisson  and  the  service  times  were  exponential  the  applicability 
of  the  analytic  approach  would  be  greatly  broadened  by  using 
M/M/N  models.  This  program  was  written  in  Fortran  IV  language  and 
resides  on  the  DLC-10  computer  system  of  'I'SC  under  the  file  name 
CHISQ.F4.  The  program  listing  is  included  in  the  appendix. 


2.3.2  Postprocessing  Program  K0LM0G.F4 

This  program  was  written  to  compare  the  model-predicted  waiting 
time  distribution  with  the  field-observed  data.  The  computed  and 
the  observed  waiting  time  distributions  are  the  inputs  to  this 
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program.  The  Kolmogroff  statistic*  is  computed  from  the  input 
data.  Based  on  this  statistic  a decision  is  made  whether  to  accept 
or  reject  the  hypothesis  that  the  predicted  waiting  time  distribu- 
tion can  be  used  to  represent  the  field  observed  waiting  time  dis- 
tributions. The  outcome  of  this  test  determines  whether  or  not  the 
queuing  model  can  be  used  to  represent  airport  landside  facilities. 

As  before,  the  program  was  written  in  Fortran  IV  and  resides 
on  the  DEC-10  computer  system  at  TSC  under  the  name  of  K0LM0G.F4. 
The  program  listing  is  also  included  in  the  appendix. 


K.V.  Bury,  Statistical  Models  in  Applied  Science,  John  Wiley  and 
Sons,  New  York,  1975. 


3.  MODEL  VALIDATION 


3.1  DEFICIENCIES  IN  THE  FIELD-OBSERVED  DATA 

Data  deficiencies  are  discussed  here  for  two  reasons.  First, 
because  of  these  imperfect  data  a judgement  of  the  validation 
results  should  be  qualified  accordingly.  Second,  if  more  data  are 
to  be  obtained  for  further  analytic  model  validation,  those 
deficiencies  are  the  pitfalls  to  be  avoided.  It  may  be  noted  that 
the  field  data  were  not  taken  specifically  for  the  purpose  of 
analytic  model  validation  but  rather  for  a simulation  validation. 
This  would  probably  account  for  some  of  the  reasons  for  a deficient 
data  set.  The  major  problems  of  the  data,  in  addition  to  the  ex- 
pected nonstat ionarity , are  listed  as  follows: 

a.  The  majority  of  arrival  data  was  recorded  as  number  of 
arrivals  per  5 minutes.  That  means  that  inter-arrival  times  can 

be  erroneous  by  as  much  as  5 minutes.  However,  because  the  service 
times  are  on  the  order  of  seconds,  there  is  an  incompatibility  in 
the  accuracy  of  the  two  sets  of  data.  A 5-minute  inaccuracy  in 
data  is  too  coarse  to  draw  definitive  conclusions. 

Waiting  times  were  grouped  in  5-second  intervals  to  form 
distributions.  Again,  compared  to  the  expected  service  time  this 
is  also  too  coarse  for  detailed  validation.  It  is  suggested  that 
the  arrival  times  and  the  waiting  times  for  each  customer  should 
be  recorded  and  aggregated  in  1-second  intervals. 

b.  A variable  number  of  service  channels  were  used  during 
the  data-taking  period.  The  number  of  servers  changes  as  the 
demand  slackened  or  increased. 

c.  The  data  are  incomplete.  For  all  the  facilities,  the 
queue  length  data  are  missing. 

d.  Observed  data  were  arbitrarily  sampled.  This  may  or  may 
not  give  a reasonable  representation  of  the  population  on  which 
the  data  were  taken. 


In  spite  of  these  deficiencies,  three  facilities  were  selected 
for  model  validation.  They  are  the  UAL  express  check-in,  the 
security  check,  and  the  Branniff  gate. 


3.2  TEST  PROCEDURE 

The  field  observed  arrival  data  and  the  service  time  data  were 
preprocessed  to  conform  with  the  input  requirements  of  the  G/G/l 
model.  Based  on  the  preprocessed  data  the  waiting  time  distribu- 
tion and  the  queue  length  distribution  are  computed  by  the  model. 
The  predicted  waiting  time  distributions  and  their  expected  values 
are  compared  with  the  field-observed  counterparts.  The  comparison 
of  the  mean  is  done  by  the  t-test  and  the  Kolmogroff  test  is  used 
for  the  distribution. 

3.2.1  T-Test  for  Means 

The  mean  value  of  the  predicted  waiting  time  is  compared  with 
the  observed  mean  waiting  time  using  the  t-test.  The  null  hypoth- 
esis for  this  test  is  that  the  population  mean  of  the  observed 
waiting  times  is  equal  to  the  predicted  value  on  the  basis  of  a 
random  sample.  The  t-statistic  is  defined  as: 


s//n 

where 

x = sample  mean 
a = predicted  value  of  mean 
s = standard  deviation  of  sample 
n = number  of  sample  points 

I'lte  absolute  value  of  the  t-statistic  is  compared  with  a critical 

value  t /,  for  a given  significance  level  a.  The  null 

hypothesis  will  not  be  rejected  (i.e.,  it  will  be  accepted  that  the 

mean  waiting  time  equals  a)  if  the  absolute  value  of  t-statistic 

is  less  than  t ..  The  numerical  values  of  this  critical 

a/2,  n- 1 
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statistic  may  be  found  from  standard  statistical  tables.*  The 
significance  level  a is  usually  set  at  0.05. 


3.2.2  Kolmogrof f Test  for  Distributions 

For  a given  sample  x of  sire  n,  such  that 

1 - 2 - n 

a cumulative  distribution  function,  G^fx),  is  defined  by  the 
fol lowing : 

Cn(x)  = 0,  x < Xj 


= r/n,  xr  < x < x^ 


1 ’ x - xn  ' 


If  the  sample  x comes  from  the  completely  specified  postulate, 
F (x)  , then 


Pr  i G fx)  - F(x)  = 0 


for  all  x . 


Based  on  this  fact,  a test  of  fit  for  F(x)  is  constructed  by 
defining  the  Kolmogroff  statistic. 


D = sup  G (x)  - F(x) 
x 


If  the  value  Dn  is  larger  than  a critical  value  (as  a function  of 
n)  (dn)c,  the  null  hypothesis  that  the  sample  x comes  from  the 
postulate  F (x)  will  be  rejected.  For  significance  level  a = 0.05 

4* 

the  critical  value  is 


rSee,  for  example,  E.I..  Crow,  F . A . Davis,  and  M.W,  Maxfield, 
Statistics  Manual,  Dover  Publications,  New  York,  1955. 

^Bury,  Statistics  Models,  op.  cit. 


A 


(d)r  = , for  a = 0.05  . 

n c /n 

The  observed  waiting  time  distributions  are  tested  by  the 
Kolmogroff  test. 

3.3  FACILITIES  TESTED 

3.3.1  Security  Check 

3.3. 1.1  Observed  Data  - Data  were  taken  on  January  23,  1976, 
between  1500  and  2000  hours.  The  arrivals  at  this  station  were 
aggregated  over  5-minute  intervals,  which  resulted  in  59  data 
points.  From  these,  the  histogram  shown  as  Figure  1 was  produced. 

A chi-square  test  of  the  data  showed  no  agreement  with  a Poisson 
process,  and  thus  the  G/G/l  model  was  applied.  Two  channels  were 
generally  in  service;  however,  at  infrequent  intervals  only  one 
channel  was  operating.  The  service  time  per  channel  was  operating. 
The  service  time  per  channel  was  assumed  to  be  a constant  8 
seconds  per  customer.  Waiting  times  were  recorded  for  726  arrivals 
sampled  from  a total  of  3,511  actual  patrons.  The  frequency  dis- 
tribution of  waiting  times  is  shown  in  Table  1. 

3. 3.  1.2  Test  Results  - In  order  to  get  better  resolutions,  the 
arrival  rate  was  converted  from  the  number  of  arrivals  per  5 
minutes  to  number  of  arrivals  per  4 seconds.  This  is  done  by 
dividing  the  observed  number  of  arrivals  (per  5 minutes)  by  75. 

It  is  true  that  the  resulting  arrivals  may  not  represent  the  actual 
arrivals  if  better  resolution  were  used  when  the  data  were  taken. 
However,  this  is  probably  the  best  that  can  be  done  with  the  exist- 
ing data. 

In  the  majority  of  times,  two  channels  were  in  operation, 
each  with  a constant  8 seconds  service  time  (estimated).  The 
service  time  used  in  the  analytic  model  was  assumed  to  be  4 seconds 
for  97  percent  of  time,  12  seconds  for  2 percent,  and  20  seconds  for 
I percent  of  time.  It  is  believed  that  this  service  time  distribu- 
tion used  in  the  analytic  model  validation  is  more  realistic  than 
the  estimated  constant  service  time. 
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IGURE  1.  OBSERVED  ARRIVALS  AT  SECURITY  CHECK 


TABLE  L.  OBSERVED  WAITING  TIME  (IN  SECONDS)  AT  SECURITY  CHECK 
FRIDAY,  JANUARY  23,  1976,  1500-2000  HOURS 


The  comparison  of  mean  and  distribution  are  done  separately 
us  follows: 

a.  Test  for  Mean  Wating  Time:  The  observed  11.63  seconds 

mean  waiting  time  at  security  check  is  compared  with  the  model- 
predicted  mean  waiting  time  of  10.41  seconds  by  the  t-test.  The 
t-statistic  is  computed  as 


The  critical  value  is 


11.63 


10.41 


23.05/  JITS 


= 1.425 


ta/2,n-l  = t0. 025,725 


1.960  . 


Since  t is  less  than  ta/2  n-1’  t*iere  i-s  n0  rcason  to  reject  the 
hypothesis  that  the  G/G/l  model  can  be  used  to  predict  the  mean 
waiting  time  for  the  security  check  tested. 

b.  Test  for  Distribution:  The  observed  and  the  computed 

waiting  time  distributions  are  compared  by  Kolmogroff  test.  The 
Kolmogroif  statistic  is  computed  by  the  program  K0LM0G.F4  as 
follows : 


Dn  = 0.08099,  with  n = 726 


The  critical  value  is 


(d  ) 
v n c 


1.36  1.36 


= 0.055,  for  n = 0.05 


Since  D is  greater  than  (dn)c>  it  is  concluded  that  the  G/G/l 
model  cannot  be  used  to  predict  the  waiting  time  distribution  at 
the  security  check  tested. 


The  field  observed  and  the  computed  waiting  time  distributions 
are  plotted  in  Figure  2.  It  is  noticed  that  for  waiting  time  less 
than  10  seconds  the  field  observed  curve  is  flat.  This  flat  por- 
tion is  the  cause  of  the  high  value  of  Kolmogroff  statistic  Dn. 
Checking  with  the  observed  raw  data  it  is  found  that  there  are 
many  (505  out  of  726)  who  went  through  security  check  without  any 
wait.  All  of  the  rest  waited  10  seconds  or  more.  However,  there 
were  no  persons  recorded  as  waiting  less  than  10  seconds  if  they 
were  waiting  in  line.  This  most  likely  indicates  that  the  person 
who  took  the  data  arbitrarily  recorded  anything  less  than  10 
seconds  as  0 seconds.  This  casts  some  doubt  on  the  validity  of 
this  portion  of  the  field  data,  which  therefore  was  omitted  in 
comparison.  The  observed  and  the  computed  results  are  summarized 
in  Table  2. 

3.3.2  UAL  hxpress  Check-In 

3.3.2.  1 Observed  Data  - Data  were  taken  on  January  23  , 1976, 
between  1500  and  2000  hours.  The  arrivals  were  taken  in  5-minute 
intervals.  The  histogram  of  the  arrivals  is  shown  in  Figure  3. 

The  arrivals  do  not  agree  with  a Poisson  process. 

A total  of  58  arrival  intervals  were  observed.  The  number 
of  channels  in  service  varied  between  one  and  two.  The  service 
time  per  channel  was  in  cumulative  distribution  form.  (See 
Figure  4.)  There  were  192  arrivals  sampled  from  a total  of  337 
actual  passengers.  (See  Table  3 for  details.) 

3. 3. 2. 2 Test  Results  - As  before,  the  arrival  rate  was  converted 
from  number  of  arrivals  per  5 minutes  to  number  of  arrivals  per 
10  seconds  for  better  resolution.  The  service  time  distribution 
was  deduced  from  the  observed  cumulative  distribution  plot  (Figure 
4 ) . 

The  effect  of  the  variable  number  of  channels  in  operation  was 
compensated  by  multiplying  the  arrival  rate  by  a factor  correspond- 
ing to  the  average  number  of  servers  in  operation.  The  comparison 
of  mean  and  distribution  are  done  separately  as  follows: 


1 5 


CUMULATIVE  PROBABILITY  OF  WAITING  TIME 


FIGURE  2.  SECURITY  CHECK  WAITING  TIME  DISTRIBUTION 
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TABLE  2.  SECURITY  CHECK  COMPARISON 
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EXPRESS  CHECK-IN  SERVICE  TIME  DISTRIBUTION 


TABLE  3.  OBSERVED  WAITING  TIME  (IN  SECONDS)  AT  UAL  EXPRESS 
CHECK-IN,  FRIDAY,  JANUARY  23,  1976,  1500-2000  HOURS 


ML  I. 
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a.  Test  for  Mean  Waiting  Time:  The  observed  59.75  seconds 

mean  waiting  time  at  UAL  Express  check-in  is  compared  with  the 
model -predicted  mean  waiting  time  of  52.68  seconds  by  the  t-test. 

The  t-statistic  is  computed  as: 

t = *— ± 
s//n 

59.75  - 52.68 
68.15/  /T5T 

= 1.437. 

The  critical  value  is 

ta/2,n-l  = t0. 025, 191  = 11960 

Since  t is  less  than  t(X/,2  n_i>  there  is  no  reason  to  reject  that  the 
H/G/l  model  can  be  used  to  predict  the  mean  waiting  time  for  the 
UAL  express  check-in  tested. 

b.  Test  for  Distribution:  The  observed  and  the  computed 

waiting  time  distributions  are  compared  by  Kolmogroff  test.  The 
Kolmogroff  statistic  is  computed  by  the  program  K0LM0G.F4  as 
follows : 


D = 0.11894,  with  n = 192 
n 

The  critical  value  is 


<dn>c 


1 . 36  = 1.36 
✓n  /T92 


= 0.098,  for  a = 0.05 


Since  Dn  is  greater  than  (dn)  , it  is  concluded  that  the  fi/G/1 
model  cannot  he  used  to  predict  the  waiting  time  distribution  at 
the  UAL  express  check-in  tested. 

The  field  observed  and  the  computed  waiting  time  distribu- 
tions are  plotted  in  Figure  5.  Note  that  if  the  waiting  times 
which  are  less  than  30  seconds  are  treated  the  same,  the  predicted 
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CUMULATIVE  PROBABILITY  OE  WAITING  TIME 


waiting  time  distribution  is  a rather  close  fit  to  the  observed 
waiting  times. 

The  observed  and  the  computed  results  are  summarized  in  Table 
4 for  better  comprehension. 

3.3.3  B r a n i if  Gate 

3.3.3. 1 Observed  Data  - Data  were  taken  on  January  25,  1976, 
between  1000  and  1500  hours.  Forty-four  arrivals  were  observed 
and  recorded  to  the  nearest  minute.  The  histogram  of  these  is 
shown  in  Figure  6.  A chi-square  test  showed  a lack  of  agreement 
with  a Poisson  process.  The  service  time  per  channel  is  shown 
as  Figure  7 in  cumulative  distribution  form.  In  this  case,  only 
one  channel  was  in  service.  Braniff  gate  waiting  time  data  were 
combined  with  all  other  gate  waiting  times  to  form  a single  dis- 
tribution shown  in  Table  5.  There  were  368  customers  sampled 
from  a total  of  1,498.  This  combined  dsitribution  was  used  for 
comparison  with  the  G/G/l  model. 

3. 3. 3. 2 Test  Results  - For  better  resolution,  as  before,  the 
arrival  rate  was  converted  from  numbers  of  arrivals  per  minute  into 
number  of  arrivals  per  10  seconds.  The  service  time  distribution 
was  deduced  from  the  observed  cumulative  distribution  plot  (Figure 
7).  The  mean  waiting  time  predicted  by  the  G/G/l  model  at  Braniff 
gate  was  33.85  seconds,  which  is  significantly  lower  than  the 
observed  mean  waiting  time  at  the  gate  (68.96  seconds).  However, 
it  is  emphasized  that  the  field  observed  data  were  not  observed  at 
Braniff  gate  alone.  They  were  accumulated  for  all  gates  in  service 
during  the  observation  period.  Since  the  Braniff  gate  was  relative- 
ly lightly  used  and  the  service  time  was  generally  lower  than  the 
rest  of  the  gates.  The  predicted  lower  waiting  time  at  Braniff 
gate  was  consistent.  Because  of  the  lack  of  field  data,  the 
statistical  tests  were  not  done.  However,  the  observed  (for  all 
gates)  and  the  computed  (for  Braniff  gate)  results  are  summarized 

in  Table  6. 


'ABLE  4.  UAL  EXPRESS  CHECK-IN  COMPARISON 
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NUMBER  OF  ARRIVALS  PER  MINUTE 


FI  CURE  6.  OBSERVED  ARRIVALS  AT  BRAN  IFF  GATE 


25 


TABLE  5.  OBSERVED  WAITING  TIME  (IN  SECONDS)  FOR  ALL  GATES, 


99%  waited  less  than  455 


TABLE  6.  BRANIFF  GATE  COMPARISON 


4.  CONCLUSION 


As  mentioned  before,  the  field  data  were  not  collected 
specifically  for  the  analytic  model  validation  and  thus  were 
not  ideal  for  the  purpose  of  this  study.  Instead  of  taking  the 
number  of  arrivals  per  5 minutes  and  taking  the  number  of  delays 
in  a time  bracket,  the  data  for  actual  delay  in  seconds  for  each 
passenger  would  be  more  suitable  for  analytic  model  validation. 

Because  of  the  single  channel  limitations,  there  are  many 
facilities  in  the  airport  landside  which  may  not  be  represented  by 
the  G/G/l  model  (e.g.,  ticket  counter,  baggage  claims,  curb  sides). 
However,  despite  the  less  than  desirable  field  data,  for  all  three 
facilities  tested  the  G/G/l  model  adequately  predicted  mean  values 
of  the  waiting  time.  Moreover,  the  tested  facilities  include 
those  whose  numbers  of  servers  were  changed  to  accommodate  demand 
changes.  It  seems  to  suggest  that  as  far  as  the  mean  values  are 
concerned,  the  G/G/l  model  can  be  used  for  waiting  time  predictions. 
The  predictions  of  waiting  time  distributions,  however,  were  not 
good.  The  model  did  not  do  well  enough  to  pass  the  Kolmogroff 
test,  possibly  because  of  the  reasons  described  above. 
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APPENDIX 
PROGRAM  LISTINGS 

A . 1 YNEVTS.F4  - G/G/l  MODEL 


C T *- 1 5>  PROGRAM  .AS  OFvElCPEO  HY  PARCEL  F.  NEuTS  - DEF ART  - 
C mFNT  OF  STATISTICS  - pOrLLE  UNTvErSITY  - wF  ST  LAFAYETTE- 
C INDIANA.  OCTOBER  1«71. 

c 

C Thic  program  POMFcTFS  T PE  TIME  DEPENDENT  FFATURES  OF  A 
C SINGLE  SE*VFP  OISCPETt  TINE  uLFLE  WITH  A F I NI 1 E WAITING- 
C POOP. 

C to  I T P TPt  FPFSFn.T  DIMENSION  STATEMENTS  A w A l T I NGPOOM  OF 
C SI7E  LP  TO  CNF  PLNDPFl  pay  Pt  MUOIEU.  THE  DENSITY  OF  T PE 
C SFPVlCF  T I E C AN  PE  f.  0 N C E N TPATcC  ON  uP  TO  THIRTY  POINTS. 

C 

C T P£  FOLl  OUTFIT  CF  T p i R PrOGrAm  INCLUDES  ThE  FOLLOWlNGv 
C 

C 1.  TPE  m F a n UlcUt  L E N 0 T P AI  T I t N. 

C ?.  T nf  l,lsTWIcUT  1CN  OE  TPE  CuFut  LENG1H  AT  TIME  N. 

C I pE  mean  «(I  T INCT  Imp.  at  uh  N. 

C A.  1 HE  DISTRIBUTION  FE  1 P c to  A I T 1 NGT  IMF  AT  TIME  N. 

C P.  THE  jC1kT  DENSITY  lE  TPl  CLFLE  LENGTH  Ac U THE  RESIDUAL 
f.  FtPVlCE  TPf  AT  TIME  N. 

C 

C At  L THESE  APE  rOMP'.TE  u EC"  N U-  TO  A SPECIFIED  VAHJE  NNN  . 

C 

C PY  USE  (,F  TPE  V«'P  IOLc  UPlIuNS.  LISTED  PELO-  ♦ SOME  OF  THESE 
C FEATURES  m A Y j(-  L E L E T t L EhOM  T w£  to  p I T T E N OUTPUT. 

c 

C THE  T *- 1 f . P E T I f ; a | LEVELpPMENl  uf  1 HE  LISCPtTF  TIME  QUEUE 
C WIlP  Af  UNHOUM/EL  OLE  lE  LENGTH  may  tt  FOUNp  INv 
( 

C * ST  E LL"  C.  UAE'EPmOS  AND  mAPCFi  F.  NEUTc 
C * A STM-LC  S E P V E P OLELE  IN  D i Q F P E T E TIME  * 

C * CAPlt^S  Dl  CEN.lPE"  Pt  "ECHExC^t  CP  t R A T 1 ON ‘ 'E  LL E — ' 1 V 7 1 . 
f 

C TPE  F (.  E | 0 » T \ G I A CUM  PAN  i UN  P a P E.  R TO  TPE.  ‘■•RESENT  pRCCRAMv 
C 

C * mArCE  L F . N Ft,T  s 

C » TPE  ^ T .m  i - 1 E <--EPVt'P  f L E L t IN  FtsCpeTe  TIME  - N UME  P 1 L AL 
t « ANAl.  Y s I s « 

C « PLPDUt  i"  I m F F Ur  A P p St PIES  - PFPT.  OF  SIATTSTICS. 

C * PURDUE  UN.  I V c I-  S l 1 Y - toE  s I LAP  Y 1 1 IE  - IN  - A7907 
C 

1. 1- c NS  TON  r ( 1 o'l ) .r  ( if; ) «PH  ( l Oil  «3c  ) « X ( llK' ) . Y ( 100  . JO) 

I 1 ME  NS  1 ON  7 ( lo  n ) . t.  ( l 0 0 ) , n T ( j 0 U C ) 

INTE  r-fcp  (J^TdC) 
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o o r->  r~>  o 


ttlHHIv  T F ( C H I I C N c « #<l<l  u 

IN  CRCEp  TP  biCITt  OUT  ll-t  JC[M  OtNSllY  Ob  1 Hfc  OUt lit  LENGTH 
AM)  ThF  RtSlDI  al  SERVICE  Ui«t,  StT  CPTT1>=1  - OTHERWISE 
SET  OPT (])=» 

READ  (H.QRc  ) C R T (It 

IN  ORDER  TO  went  CL  I IRE  «_  li>T'lrtOT  ION  CF  THE  -J  A I T 1 Nb- 
TImE.  SET  OPT  ( P ) = 1 - LlhtFuIjf  SKI  CPT(E)=e 

PE  AC (S.PRR)  CFT  (P) 

THE  US  Eh  may  '-TSH  TC  COMRLIt  Art  WHITE  THE  C1STRIPUTION  OF 
THE  WAlTINOTln  CM  Y wT  T I f-.  fc  POINTS  WHICH  ARE  A MULTIPLE  OF 
A CONSTANT  M.'.  THIS  To  SAVE  UK  RPUCEbSlNb  TIME  AND  ON  Tnt 
N U M P t R OF  lli-OH  c F C U I P L I . IN  THIS.  CASt  THE  IDENTIFIER 
OPT  (3)  SRUUlL  OF  SET  euLAL  To  nNE  AM}  THE  NUMb  F. R NR  ShOLLU  bK 
C-IVFN.  CP1(3)  AM  NP  AHt  To  hF  GIVEN  IN  AN  II.IF  FORMAT. 

REAL  ( b • PPM  CPT  (3)  .NR 
XNR=NR 

KTE  ST  I =OPT ( 1 ) 

KT  EST?  = PPT  {?) 

K T E ST  3 = 0H  T ( 3 ) 

IF (PTfST3.F0.1)  k I f SI E=0 
IF (KTEST3.F0.0)  K|tSlA=0 


eiHtt « Y h F 
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LI  IS  ThF  SI2F  OF  THE  V;  A IT  iNuRf  CM . Ll  IS  LT  . 101 
HEALTH. 10^1)  Ll 

L 2 IS  I HF  M l M D F R CE  PCIMS  ON  MICE  THF  DENSITY  OK  T Ht 
SERVICE  TIME  TS  CONCENTRATED.  | c IS  Ll.  31 

REACTS.] 00  I)  L? 

K IS  THE  MAXIMUM  N U m H t R CF  ARP  T V AL  S PEp  UNI  1 CF  TIME. 

K SHOULD  hE  AT  LEAST  CNE  AND  STPICTlT  LEsS  Than  Ll. 

REACTS.  10  0 1)  K 

RTj)  IS  THE  PPOHAHIlITY  Teat  a SERVICE  TIME  LASTS  FOR  J 
UNITS  OF  TIME. 

ONE  SHOULD  V E r-  I E Y PEFcpERANE  ThAT  THE  SUM  e>  ( 1 ) ♦ . . . R ( L £ ) 
IS  FQUAL  TO  One. 

REACTS. 10O?)  (R  ( J)  , j=l  .Li?) 


P(J)  IS  THF  PPOPAblLITY  THAT  j CUSTOMERS  JOIN  THE  QUEUE 

DURING  A UNIT  OF  T I mF . 

the  index  j pins  from  one  to  k. 
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READ(5*100?)  (P(j)*j=1.K) 

C 

C PO  IS  ThE  PHCFABlLlTY  THAI  NO  CUSTOMERS  ARRIVE  DURING  A 
C UNIT  OF  Tlf'F. 

C ONE  SHOULD  VFCIFY  RFFdREHANC  That  PO  ♦ P(l)  ♦ . . • >P  < K > 

C IS  EQUAL  TO  Ok  F . 

c 

REAL  (Si  1 0^:1)  PO 
C 

C Id  IS  the  initial  OLFlE  LENGTH. 

C JO  IS  THE  initial  PFSIOLAl  SEM/ICE  TINE. 

C IF  10=0*  THFN  jo=c  and  conversely. 

C 10  SHOULD  NOT  FXCEEC  LI. 

C JO  SHOULD  NOT  EXCEED  l2. 

c 

RFAD(S*1001>  10 
HEAD (b,  10  01 ) JO 
C 

C NNN  IS  THE  MAXIMUM  TINE  HCINT  FCR  WFICH  THE  QUEUE 
C FEATURES  ARE  HORFLTFD.  NNN  SHOULD  FE  AT  LEAST  ONE  AND 
C at  HOST  9999.  NOTE  HOwEvER  THAT  THE  PROCESSING  TIME 

c and  the  number  of  lines  of  output  crow  proportionately  to 
c the  value  of  nnn . 
c 

HE  AD ( S t 1 0 OR ) NNN 
IF (lO.EQ.O.AND.oO.EG.O)GUTn  2001 
PP<10*J0)  = ] . 

GOTO  2002 

2001  P00=1. 

2002  N = 0 

XI  l=PO 

DO  21  1=1 .H 
X I = I 

X 1 1=X 1 1 *R  ( I) 

X 1 = X 1 ♦ X 1 »p  ( I ) 

21  CONTINUF 

DO  22  J= 1 «L2 
Xj  = j 

X11=X11*R(J) 
x2=x2*xj*r ( j» 

22  CONTINUE 
X11=X11-?. 

XII  = AtiS  (XI  1) 

WRITE (3* 10 1A) 

WRITE  (3.1022 )K*L2*Ll 
N1  =0 

wRITE(3.1007>M.P0.(j*P(J>,j  = 1»K) 

WRITE (3. 1012) 

WRITE  (3*1010)  (J*R(J)  t'w  = l*L?) 
wR  HE  (3*1012) 

WRITE (3*1008) 10* JO 
WRITE  (3.1012) 

WRITE  (3*  1009) NNN 
WRITE (3.1012) 

WRITE (3*1013>X1*X2 
WRITE  (3*1006) 
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c p 1 a C-  N 0 S T i c r **««* 

c 

IF (Ll.LE.K)  GCTC  cOOS 

IF < IO.FO.O.ANC.jO.NE.Q.OR. IO.NE .0. AND. jO.EQ.OIGOTO  2 0 05 
IF  (Xll.C-T..  000001  ) GOTO  20o5 
IF (KTEST3.EC.1  .ANL.NH.LT .2)  GC1C  2005 
WRITE (A, 3200) 

3200  FORMAT  </»lX.*THF  ►-ROC-RAM  (-  a S RCN# ) 


LI  1 =L 1 - 1 
L2 1 =L2- 1 
L22=L2 ♦ 1 
M I =L 1 *L  ? 

K 1 =K ♦ 1 

UO  13  1 = 1. LI 
Y ( 1 «L2>  =0.0 
13  W ( I ) =1 .-Pr 

IF  <K .EO. 1 1 GOT  C 2003 
DC  eO  1 = 2. K 
W ( I ) =W ( T - 1 ) -P ( I - 1 ) 

20  COM  I NOE 

AT  This  STAGE  THE  INPuT  DATA  (- A V t SEEN  REAn  IK,  THE 
pp-ARRAY  h-At  PEEK  INITIALIZE!!  AND  THE  IKPLiT  DATA  HAVE 
PEEK  WRITTEN  GUT  AKC  ^UbJECTET  TO  SOME  RUDIMENTARY 
DIAGNOSTIC  TESTS.  THE  NEX)  LINE  STARTS  THE  MAIN  LOOP 
WHICH  IS  REPEATED  KKK  TIMES. 

otMMMt  the  MAIN  LOoP■,M,',>^‘,► 

2 0 03  Mi  = N ♦ 1 

IF  (K.uT  .M  Ki ) STOP 
C 

C THIS  PORTION  OF  IHP  PROGRAM  COMPUTES  the  new  PP-ARRAY. 

C PP(I.j)  IS  T HP  PRC P ARIL  ITY  that  Al  1 HE  TiMF  CCNSIDEREU 
C THFRE  A h E I Cl  STCMF-R  IN  the  SYSTEM  ANU  The  RESIDUAL 
C SERVICE  TIMt  OE  IRE  CUSTOMER  PpING  SERVED  IS  J. 

C THIS  IS  FOD  I MFTwEFN  ONE  Anu  I 1,  FCR  J bETWEEN  ONE  AND 
C L 2 . TRt  IDE  NT  TF  ItR  ROu  CGMAINc  Th£  PRObARiLITY  THAT  THE 
C OuEOE  Is  E^RTv. 

C 

WR II E ( 3, 1 GOO) 

WRITE ( 3.  IP]  7)N 
1»R  IT E ( 3 . 1 A I 2 ) 

XGG=PUO*po  (1,1) 

DC  1 1 = 1.* 

x ( 1 ) =R  (T)  *R00*P'j*HP  ( I ♦!  , 1) 

11=1-1 

DO  3 NU=t.I 

X ( 1 ) =x  ( I ) ( 1 ♦ 1 -N'L  ) «RR  INU.  1 ) 

3 CON  I I NI  T 
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u 1)0  « JM.I21 

Y ( I , J) =PC«PF ( l , j*  1 ) 

IF (l.tO.l)GClC  51 
DO  5 NV  = 1 . 1 1 

Y (I  , J ) = Y (T  » J ) ♦D  (T-NV)flPP(N\/.J4l) 

5 CONTINUE 

51  COM  I NUF 
? COMINUF 
1 COMINUF 

DC  6 I=K 1 ,L  11 
N01=I-K 
UU2  = NU 1 ♦ 1 
I 1 = 1-1 

X ( 1 ) =fO*Pf  (1*1.1) 

DO  V NU- M H . 1 

X(I|=A(I)*F(l*l-M.)»uMM).U 
7 COMINUF 

DO  S.  j=l,|.?l 
Y ( I . J) =PO»PP  ( I . j*  1) 

(JO  H M'  = NI  1,11 

Y(1,J)=Y(T,J)*F(T-NV)«PP(N\/.J«1) 
f COMINUF 
9 CONTINUE 
5 CCNUM'E 
X ( I.  1 ) = 0 . 0 
00  1 0 * 0= 1 ,K 

X(L1I=<(L1)*MM.)«P';(L1  + 1-»L,1) 

I u COM  I NUF 

I/O  11  J=1.L2J 

Y (L  1 » o) =Pp (L I . J*1  ) 

no  \?  M'=i  .k 

Y (l.  1 »o)  =Y  (l  1 . w)  ♦I'lNU  ) *PP  <L  ) -M'liv*  1 ) 

1e?  CCNTIM'F 

II  CONTINUE 
POO  = PU*»G'' 

c 

C 7(1)  CONTMM  E IK-.T  Trt  UENM'lY  ANC  NF  XT  T >-t  CISTWIhUT  ION 
C rf  T (.Ur.OF  it  NOTE-  A i f*-t  TIM  FC1M  CONSfOF>tO. 

C XVN  CCOlAlM  TEE  I/*  AN  .UtLE.  LFU(:TE  AT  T Ft  TImE  POINT 
C CCNSIOt fE d. 
c 

X M i = 1 .-PC" 

Of.  1 A 1 r 1 ,1  1 

/ ( 1 ) =.).u 

('(.  15  .1=1.  L 7 

P';  ( ! • J)  =Y  (1  <„  I <’X  ( 1 ) 

/ < T ) =/  ( I)  te  a < I .-j  ) 
lb  CO'-.TINO*- 
1 <♦  ( ri  I l i">jh 

/ ( 1 ) = / ( 1 ) ♦ 1 C i 
>w-t'..sXN"f  ♦ 1 ,-l  ( 1 ) 

00  I F l:E.l  1 
l ( 1) =/ I 1)  .7  ( 1-  1 ) 
xm.  =xhf  ♦ 1 ( | ) 

1 a COM  INOE 
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IF  ( k T F S T # E u • o ) CM  ( ; 2 0* jj5 
XN  = N 

l ) J = « N1  Ij  i • l X v . X N p ) 

I F ( 0 i . L T . ■'  . A ) K T E s T a = 1 
IF  ( 1. 1 . r-T  . i . a ) MF  3 U = C 
IF  (i'TrcT‘..c,;,l)  Ot.  I 0 >Ulo 
J F ( F 1 2 c T <i  . ,r  C • i; ) e i<  u A 

2U0o  CCMIM.'F 

THIS  H.JkIIon  f'f  7hF  Pr<u0f-AN  CC'-PUTf.S  1 i-h  C I ST  F.  IHUT  1 UN  uF 
T H t vlelU^L  wAITINGTImF  a 1 TIN*-"  N • 1 6 F ALGORITHM  I1-  AN 
ANALOtUt  OF  hrKf,tH*c  k t T hCl<  EC  p T *-E  EvAIUATION  OF  QH l- 1 - 
xakY  F-OLVNOV I MS.  hlT  ADAPTED  ‘-•ckF.  10  CONVOLUTION  PkOuDCTS. 

IFtKle.STr’.FG.u)  F-  M U 2 0 0 6 
2010  COM  IM  F 

cr.  32  j=  1.1.2 
wt  t j>  =pf  (i  i > 

32  CONTINUE 

LG  3 3 J=L'2.I-1 
wT  <J)=0." 

33  CONTINUE 
MM  = 1 
MN3  = 2 

M N a = L 2 

UC  3*4  JX  = 1 .Ll  1 
JX  1 =L  t - jx 
VN2=MNA 
MN*4=MK?*L? 
mNS=MNA.? 

DO  33  J=2 . N NA 
JP  = N N3~ J 
wT  ( jP ) = 0 • o 
MNt =maxo ( 1 , jP-vn2  ) 

MN7  = mr  0(1  7,vF- i ) 

DC  33  NU  = >,N6.MN7 

e/T  ( j k ) = I*  T ( J P ) ♦ P ( N o ) <l  w t ( J K ~ * L ) 

33  COM  I NDF 
WT  ( 1 ) =n.O 
DC  37  JJ= 1 » L 2 
WT  (jj) =WT  ( J J ) ♦ P P (jXl  iwj) 

37  CONTINUE 
3a  CCNTINDF 

aT  ( 1 ) = POO*M  < 1 > 

Zmn  = 2 ,-POO-vt  1 ( 1 ) 

DC  36  J = ?.M 
w T ( j ) = w T (j)  ♦ » T (j-i) 

/V\  = ?tof"+  I ,-M  ( j) 

36  CONTINUE 
2006  CCNTINUF 


WRITE  STATEMENTS  F( 

WP  I T E ( 3 « 1 0 1 6 ) N • X M|\ 
WRITE  (3. 1 '12) 


FFCUlrEO  OUTPUT. 
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white  n.iniP)N 
WHITE  (7.1  12) 

DC  *000  1 = 1. LI 

IF  (2(I).LF.0.9999'7)  GC  TO  *oOO 
L 11-'=  I ♦ 1 
GO  TO  *010 
*000  CONTINUE 
*0  10  CCMINUF 

«nF*ITt(.T.1011)M.cl)0»(l»2(I)»I  = l»LlP) 

WWI  It  (^.lf.l?) 

IF  <KTE5T2.F<V.0.ANL.KTEST*.rC.0>  GOTO  2007 
hk  I ( t ( '*  . 1 Or1 1 ) N » ? Mi\ 
wkl  It  r».lri2) 

GO  *lu"  1 = 1 • n 1 

IF  (WT  < 1 > .1  <■  . 0.95999)  GOTO  *100 
U 1 P = I ♦ 1 
GO  10  * 1 1 0 
*100  CONTINUE 
*110  CONTINUE 

WPITE(3,ln?0)N.M,POU»(j,WT(J)»J=l.MlP) 

2007  CCNTINUF 

IF  (KTcST  1 .EC.O)  GOTO  cOO* 
wHITt  (3.1M2) 

WHITE  0.1  0 1*5)  N 
WHlTEO.lClblFOO 
WHITE ( 3 , 1 0 1 2 ) 

00  17  1 = 1.  LI 

WH  ITF  O.  loth)  1 . (P*  < I ♦„  > » J=1  »L2> 

17  CONTINUE 
200*  UG  TO  ?00-» 

2005  WHITE  ( T.95P) 

WHITE  0.3000)  * 1 1 

3000  FCHNAT ( / . i * . »x  1 1 = F.E20.10) 

WHITE!*. 31UU)X11 

3100  FQkwAT (/.ix.xEPPGh,  *11  = r.E2C.10> 

C 

C THF  FOHHAT  ft  ATEMENTS, 

c 

*59  7 FORMAT  (11.12) 

99P  F0PmAT<*  ATTENTION*  THERF.  AkE  ENkOPS  IN  THE  INPUT*. 

•*  DATA.  Pi  EASE  ChtCK.  *) 

999  F C P n A 1 (Hi 

1000  FORMAT <*1*) 

1001  F OkM  A T ( T 3 ) 

1002  FCkMAT  (3F  7.5) 

1003  FOFwAT (F  7.*) 

100*  F CkN A T ( I*) 

10 05  FORMAT  (3X.I3.10F7,*. (EX, 1 OF 7.*)) 

loch  forma; (//i 
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1007  ECRM>AT(A  THE.  0FNS1TY  OF  ThE  N0ME1ER  OF  ARRIVALS  PER#. 

«#  UNIT  Of  TIMEa,//<2X,10<IA»F6.S>>> 

100H  FORMA1IA  THE  INITIAL  OOELF  LENGTH  IS  #.I3./#  THE#. 

«*  IM1IAL  RESIDUAL  SERVICE  TIME  IS  # » l 3 ) 

1 0 0S»  FORMAT  ( a THE  MjMRtH  CF  T IMF  POINTS  COMPUTED  IS#. 

* 1A  > 

1010  FORMA  T ( A THE  CFNS1TY  CF  1RF  SERVICE  TImES#. //2X. 

* ( 10 ( IA.FB.N1 ) ) 

1011  FORMAT  <3X.10dA,Fe.S)  ) 

1012  FORMAT  (/ > 

1013  FORMAT (2X.ATFE  m£aN  NR . uF  ARRIVALS  PER  UNIT-TIME/#, 
*F10.A,/a  THE  MEAN  SERVICE  T I ME v # , F 1 0 . A ) 

1 0 1 A FORMAT (AlA. ////A  1 HE  I RANS IE  NT  bEHAVlOR  OF  A #, 
“ADISCrME  TIME  r.OtUE  *ITH  a FINITE  V*  A I T I NGROOM  # « 

*//  » 

1 0 1 5 FORMAT  (#  the  CUtOt  IS  EMPTY  WITH  PROBABILITY#, 

*FR.b) 

1016  FORMAT  (A  AT  T I MF  N =*,IA,A  1 HE  MEAT.  OU^UE  LENGTH#. 

*#  E COALS# • F 1 C • A ) 

1017  FORMAlt#  The  COE  Oc  ChARAOTfR 1ST  ICS  AI  TIME  N = A, 

* I a ) 

1016  FORMAT (A  The  CISTkIhOTICN  hF  the  uoeue  length  a. 

* A A T T IMf  * = A . 1 H/ ) 

1 0 1 R FORMAT  (A  THt  jOIM  DENSITY  CF  THE  UOELE  LENGTH  A. 

* A AND  the  t-ESICOAL  SEhvICE  TIME  AT  TIME  N =A,lA/> 

1020  F 0 h h A T (A  THfc  DISTnlbOllON  nF  THE  RAITInC-TImE  AT  A, 
»#TlMt  f = *. 1A,// ( ja, 10 ( IA.F8.S  ) ) ) 

1021  ECh"AI(A  Tnt  MEAN  wAHINGTIME  Al  TIME  N =A.lA. 

«#  IS#  «F 1 2. A ) 

1022  E CRM A T ( a THE  LPRE*  LIMIT  CF  THE  NUMBER  OF  ARRIVALSa. 
«A  PEw  UNIT  CF  TIME  ISA. 13. /%  THE  UPPER  LIMIT  OF  ThEA. 
*A  NOMhER  CE  LMTS  OF  StRvlCE-T  IMt  PER  CUSTOMER  IS#. 

* 1 3 » / # THE  UPPER  LIMI1  1C  ThE  NUMBER  UF  CUSTOMERS#. 

«A  IN  T hf  cySTEm  IS#,  1 A / ) 

ENO 


J. 


A. 2 CHISQ.F4  - PREPROCESSING  PROGRAM 

Cl.  HEAD  IN  THE  SAMPLE  IE  APRIVAi.S  PEP  UNIT  TIME. 

C2.  SORT  THE  SAwPLt  PY  NOMbER  OF  ARRIVALS. 

C 3 . CHECK  IF  THE  SAMPLE  IS  A POlsSON  PROCESS. 

DIMENSION  A11101 ) .PI  ( 101  ) . TP < 10) 

C.  A 1 ( 1 ) =NUMHF R OF  CCCUhANCES  OF  1-1  ARRIVALS  PER  UNIT  TIME. 
C.  N = TOTAl  SAMPLE  PCIMi. 

READ(20,IOL')N 

SbM=0. 

SbM2  = (J  . 

JMAA=U 

SMALL=0.0?«N 
DC  S00  1=1.101 
SOO  Al<I)=o. 

DO  1000  1 = 1, N 
head  UO,  ICO)  1a 
IF  (1A.ME  .")f:C  TC  *001 
A1  ( 1 ) =A1  ( 1 ) +1 
GC  TO  10  0 0 
AOOO  CONTI  NOE 

DO  SOO'J  J=1.10C 
IF  ( IA.GT.jlGC  TC  SOOl 
A1(j*1)=A1  ( J ♦ 1 ) ♦ 1 
SUM  = sum ♦ j a 

SbM2=SUM?. ] A« I A 
IE ( j.GT . jv ax ) oMAXs j 
bC  IC  100° 

soon  con  iiw  f 
100C  COMIND* 

F Mt  an  = Sii«/,/i'. 
v AP  = SbM2/'  -F M E AN«E me  AN 
A A 1 |0  = VAU/EMtAN 
T Y Pp  irj.FMf  AN,  VAh.RAUC 

IE  (KA|  10.1  T.C.“.0„.hA1  IO.CT.I.cIGO  TO  PS9S 
C.  verify  lh  Tt~r  sample  is  POISSON 
JM IN= 1 
E AC  1 = 1 . 

E XPArf  XP  ( -F  ME  AN  ) 

E ME  AN  I = 1 . 

P 1 ( 1 ) =E  XPA*N 
MOO  CONlIMiE 

DO  M)iM’  I=JMJN.jMriX 

e ac  i=e  f.r  r«  i 

F Mt  uN  1 =FM=  AM  «FME  AN 
P = f MF  AN  l/p  AC  I 

6000  P1(IM>=MM 

IF  (Pi  (JMAx)  .LE  .SMuLDCO  TC  620G 
JM  IK  = jmox ♦ 1 
jv  A*  = jMAX ♦ l 
DC  TO  Ml' 

a POO  CONTINUE 
DCF  = JM  A X 
P 1 L = P 1 ( I ) 

A 1L  = A 1 M ) 


39 


r 


L 


< 


Fir.  = F 1 ( Jf  t > ) 

Mr'=“)  ( J N ' * ) 

1>C  r JO  *'  1 - r « A * 

IF  (r-u  ,r-r.  .tvALL  )<;o  TO  <310 
F iL  =F  1 1 ♦*-  1 ( I ) 

■i  1 L = A 1 L ♦ f 1 (1) 

I L = 1 

LCF  = Of.F  - 1 
F 3 1 0 COM  IK"- 

IF  (FK.f-Y  .‘••‘•■A l l.)H  To  r 
F 1 1->  = ►-  I *F  1 ( j U A * - ! ) 
ai  (-  = 4 i *-  ♦ A 1 ( Jf-  AX-1  ) 

It-  = j-tuX-  ] 

EOF  =f'UF-  l 

fcj?0  IF  (F1L.P-F  .F  lr.i,F  .SMALL  >0G  10  ^ ^ 0 

fc  30  0 CONI  lM’f 
tACO  CoM  I MjF 
i:of  =0(;f  - 1 

I L 1 = J L ♦ ) 

I »- 1 = 1 M- 1 

OF- \>.~  ( 2 U II  )*■»>/»  1 L ♦ ( A 1 t--h  1 (-)«•»</!-  1 F- 
10  O'li.i'  1-11  I • 1 *-  1 

{>  Ir  = ( 1 f ♦ ( A 1 ( ] ) --  1 ( I ) ) «■;  /F-  1 ( I ) 

F b U 0 CC  N I I M :F 

T>K  <,1  ( yr  ,OCF 

';(  | I 

S"iOO  |Yt-h  '/•'(! 

SSOF-  CCnI  IMlF 
IYki  h''* 

DC  ►•Ouc  1 = 1 1 ii 

PC  "0 10  J- 1 . H 
poio  iP<j)s(T-n  *o-i 

T YPfc  VI  u • ( IF-  <„)  » J=  1 « 1 0 ) 

T V *-fe  K’O.  t IP  ( j)  ♦ 1 ) , j=l  . i C ) 

F 0 0 0 ; > F-«-  Ff'O.  I A 1 ( IF  ( J)  ♦ 1 ) . J=1  , 1 V,  ) 

100  FCP"A7(IEI 

300  FUA^uI  I/.  ’ > ,»Fvf  t,\,  vfF.V/F  = * • 3 F 1 u . 3 ) 

*400  F0>VM  (/,lXt*CM  3<  .OCF  = f.cFIG.A) 

Ac0  FUK-A  1 <//.  1 x.*-  4KvLS/r.»Alr  FPfC/CF'SVL  FpfcO  TAKE  v>,/) 
F00  F C P-  v A 

F 0 0 FOP  *AT  ( 1 X.  ) 0 (1  * .F-  . 1 ) ) 

700  FOr-FAf  <//.l  X**TPF  s ai*P  F KN  N(1  OF  F'C  T bSCN  JF> 

END 
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POSTPROCESSING  PROGRAM 
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A. 3 


c this  mr.fcctM  (.as  levelopet  hy  li  shin  yuan,  transportation 

C SYSTF.l*"'  CPNTfcc.  Put.  CAMpR  TCGE_.  MASS..  MAY  19Tb. 

Cl.  ThI<-  PROoRAM  COmParES  TE-E  c.cStRVtt  WAITING  TIME  0 I S I E<  I «GT  I ON 

Cl.  WITH  THt  RliELtlNG  mCUEL  GENERATED  WAITING  TIME  DISTRIBUTION  USING 
Cl.  KCLMObCROFF  T EM- 
IT I Mt  NS  1 0*  OHEPC (eOC)  .GFWC (EOU) 

DO  10O0  T=1.?0C 
OhF  Ml.  | I ) a 0 . 
lifK'MlIs'. 

1000  COM  T A Uc 

4o0c  t-UHMAT  (Fb.u) 

C.  INRuTS  f °Ov  CUELEING  MODEL  : 

TYPE  MJ  0 0 

H 0 0 0 FORMAT (/, IX, *No“HER  CE  PLAsStS  PY  CuEUE INC,  MOUtL  = X./.4X) 

ACCEPT  4TOO.NO 
WR IT  P ( 1 1 .4  1 0 0 ) NO 
TYhF  Mur 

h 1 o o fQRRNT (/,)*, tTYwE  In  ThE>.  One  at  a time  »ITp  -1.  AS  THE  1ST*./T 

ACCrPT  4^>J0.  TOEru  I I > . 1 = 1 »hC  ) 

WRITE  ( 1 I . 4 ? 1 0 ) IDE  ~C.  ( 1 ) *1  = 1 .NUT 
4Elu  f uxyM  (r  7,1:) 

4100  format  ( t ;< ) 

4 E 0 u Kc»nAT (r7.C ) 

C.  FACT  Or  =GGE  1>F  I NC  UN  I I T I W t /()E  cF.  h VE  C UNIT  TIME. 

T Y P r H4  n 0 

H400  FORMAT  I/.1X.FF ACTCr  = *./.4X> 

ACCEPT  4RpU.F ACT^P 
wRl Tt ( 1 1 .4000 IFACTGr 
F A(.T0R  = 1 ./FAC  TCR 
C.  INPUTS  FRO*'  OESFRVED  CATA  ! 

J=1 

SErFw=0. 

TYPE  H 30  0 

*300  EURMAT  (/,1X.*TYhE  IN  OBc  VALLE  Y FREQ  1 PAIR  AT  a time  A./,4») 
E000  CCwllNUF 

ACCEPT  4?EO.CHSV.EktC 
WRITE  < IE.mEEO  ICPsV  ,E  Rtf. 

4EP0  FORMAT (PES.O) 

IE (Cm«V.I T.-0.01 TOPErO ( J»=SEREC 
IE  ( Omc V .L T . • C . 0 l ) GC  TO  ?1)0 
IOHSV=OHCV#E  ACTOr*  l .S 
IF (lORsv.OT.j)CO  TC  E 1 00 
SFrEC.=«F  CF.O*FREO 
CO  TC  ?orc 

E 1 no  OhE "0 ( jl = SE  REC 

J= IOrSV 

SFREO=SFRE0.FrFC 
GO  TO  tOOO 
Elio  CON  T INUF 

N = SFRE.O 

SFREd=1./SEREC 

TO  ?E  0 0 Tal.ICPSv 

OPF  R(j  • I ) =OHF  RO  ( I > *sFREO 
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22 00  CONTINUE 

C.  GET  MAX  A6C  DEVIATE  FOR  TEST  : 

DMAX=0. 

DO  2300  i=2»ioesv 

IF (ObFPO ( I ) ,LE .0.000001 >Gn  TO  2300 
IF (QFRO(I) »LE. 0. 000001) QFPt(I)*l. 

DEV=ABS (ObFRC ( I ) -OF  KG  ( I ) ) 

IF (DtV.GT ,DMAX)DwAX=CEV 
2300  CONTIMJF 

FN  = N 

SORTN^OCT  <FN) 

TYKE  SuOO ,N*SQRTN  .DMAx 

SOOO  FORMAT  (./,  IX,  ,EOK  IN .DMA*  = 1 , 1 4 , F 1 0 . 5 . F 9 . S ) 
TYPE  MOO 

MOO  FORMAT  (/,  IX, XCbSERVtC  : *) 

TYPE  ^200, (CRFKO ( I ) , 1=1 , InfcSV ) 

MOO  FORMAT  ( 5FM ,b  » / ) 

TYPE  S300 

S300  FORMAT (/, IX, ^COMPETED  : *) 

IF ( IOPSV.GT .KO)NO=I0PSV 
TYPE  RPOO, (GFRG(l) ,1=1,NC) 

EM 


2.30  Copies 
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